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Abstract. We use recent low-mass star models, which 
reproduce accurately the observed sequences of various 
globular clusters, to convert the observed luminosity func- 
tions into bolometric luminosity functions. These latter 
are shown to exhibit a similar behaviour with a rising 
slope up to Afboi ~ 9, i.e. to ~ 0.2 — 0.3 M Q , and a de- 
creasing behaviour beyond this limit. We then derive mass 
functions for globular star clusters down to the bottom 
of the main sequence. These mass functions are well de- 
scribed by a slowly rising power-law dN /dm oc to - ", with 
0.5 ^ a ^ 1.5, down to ~ 0.1 M Q , independently of the 
metallicity, suggesting a rather universal behaviour of the 
cluster initial mass functions. The effects of tidal stripping 
and mass segregation are illustrated by the overabundance 
of very low-mass stars in the outer parts of NGC6397 and 
their depletion in the central parts. This analysis confirms 
that the mass function determined near the half-mass ra- 
dius has been weakly affected by external and internal 
dynamical effects and reflects relatively closely the ini- 
tial mass function. We predict luminosity functions in the 
NICMOS filters in the stellar and in the brown dwarf do- 
mains for different mass functions and metallicities. 

We apply these calculations to the determination of the 
mass function in the Galactic halo, including the DeVau- 
couleurs spheroid and the 1/r 2 dark halo. We derive the 
slope and the normalization of the spheroid mass function 
which is well described by the afore-mentioned power-law 
function with a ~ 1.7 ± 0.2 down to 0.1 Mq, although a 
slowly decreasing mass function below ~ 0.15 Mq can not 
be excluded with the data presently available. Compari- 
son with the Hubble Deep Field star counts is consistent 
with such a mass function and excludes a significant stel- 
lar population in the dark halo. This shows that essen- 
tially all the high- velocity subdwarfs observed in the solar 
neighborhood belong to the Galactic spheroid. 

Consistent analysis with recent microlensing experi- 
ments towards the LMC shows that the spheroid and 
the dark-halo stellar +brown dwarf populations represent 
at most ~ 1% of the Galactic dark matter density. This 



clearly excludes brown dwarfs and low-mass stars as sig- 
nificant dark matter candidates. 
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1. Introduction 

The derivation of reliable stellar mass-functions (MF) 
down to the hydrogen burning limit is an essential task for 
the understanding of a wide variety of major astrophysi- 
cal problems including star formation theory, stellar evolu- 
tion, dynamical evolution of globular clusters, galactic for- 
mation and evolution. Not mentioning the long-standing 
problem of the galactic missing mass. From the stellar 
point of view, the knowledge of the low-mass end of the 
MF provides essential information about the processus of 
star formation and the physics of the interiors and atmo- 
spheres of cool and dense objects. From the galactic point 
of view, the determination of the stellar MF down to the 
hydrogen burning limit is a key issue to derive the total 
amount of matter under the form of visible, stellar objects, 
which determines the mass-to-light (m/L) ratio, and to in- 
fer the contribution of substellar objects, essentially brown 
dwarfs, to the galactic mass budget. The intrinsic faintness 
of very-low- mass stars (VLMS), and the small density of 
halo objects w.r.t. the disk population (~ 1/500, see e.g. 
Bahcall, 1986; Tyson, 1988) renders the determination of 
the luminosity function (LF) of halo stars down to the 
hydrogen burning limit a tremendously difficult task. The 
determination of globular cluster (GC) luminosity func- 
tions and mass functions is a first, easier step towards 
the determination of the Population II star initial mass 
function. Moreover, the determination of the present day 
MF (PDMF) in GCs contains potential information about 
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the dynamical history, evaporation and tidal stripping of 
the clusters, and thus about the Galactic halo history. At 
last, given the large range of mctallicities covered by all 
observed GCs, it is possible to examine the dependence of 
the MF on metallicity, and to determine whether it is a 
universal property or if it is strongly affected by the metal 
content of the initial cloud. 

Recent deep photometry observations of several GCs 
with the Hubble Space Telescope (HST) yield an accurate 
determination of the bottom of the MS of these objects. 
These results extend previous ground-based observations 
to about 4 magnitudes fainter limits and allow an accu- 
rate determination of the cluster LFs down to the bottom 
of the MS, providing stringent constraints for theoretical 
models of metal-poor VLMS. The derivation of accurate 
MFs from these observed LFs is a more delicate task. The 
MF is very sensitive to the mass-luminosity relation and 
thus the derivation of reliable MFs requires the calcula- 
tions of accurate mass-luminosity relations. Direct, obser- 
vational mass determinations of VLMS are available for 
solar metallicity objects only (Popper 1980; Henry & Mc- 
Carthy 1993). No direct mass-determination has been ob- 
tained so far for metal-depleted low-mass stars. However 
the accuracy of the stellar models can be examined by con- 
fronting the theoretical color-magnitude diagrams (CMD) 
to the observed ones. Such an extensive study has been 
conducted in a previous paper (Baraffe, Chabrier, Allard 
& Hauschildt 1997; BCAH) where these authors have de- 
rived an extended set of low-mass star evolutionary models 
(0.07 < m < 0.8 M Q ) over a large range of metal-depleted 
abundances characteristic of old disk, globular cluster and 
halo populations, namely —2.0 < [M/H] < —0.5. The ac- 
curacy of these models has been demonstrated by the re- 
markable agreement (< 0.1 mag) between the theoretical 
and the observed CMDs over the entire sequences, for all 
GCs, over the entire afore-mentioned metallicity range. 
The well known wavy behaviour of the stellar main se- 
quence, which reflects genuine intrinsic properties of stel- 
lar matter, is reproduced qualitatively and quantitatively, 
for the correct magnitudes and colors (see BCAH). As 
stressed by these authors, these models are not hampered 
by any adjustable parameter (e.g. the mixing length) and 
thus the agreement between theory and observation re- 
flects the accuracy of the physics entering the models 
(equation of state and non-grey atmosphere models). This 
assesses the reliability of these models to describe the me- 
chanical and thermal properties of metal-poor low-mass 
stars down to the brown dwarf limit, and provides for the 
first time accurate relations between the mass, the age 
and observable quantities, magnitudes and colors, for dif- 
ferent metallicities. For solar-like metallicity, the accuracy 
of the mass-magnitude relation has been assessed by the 
agreement with the observed mass-spectral class (Kirk- 
patrick & Mc Carthy, 1994; Kirkpatrick et al., 1995) and 
mass-My and mass-M^ (Henry & McCarthy, 1993) rela- 



tionships (Baraffe & Chabrier, 1996; Chabrier, Baraffe & 
Plez, 1996; Chabrier & Baraffe, 1997b). 

In this paper, we use the mass-luminosity relationships 
obtained with the afore-mentioned metal-depleted models 
to derive the mass-functions characteristic of several GCs 
observed with the HST and of halo field stars. The obser- 
vations are summarized in §2 and converted into bolomet- 
ric LFs. In §3, we derive the MFs for the different clusters 
and we examine the dependence of these MFs upon metal- 
licity, location in the cluster and Galactic position of the 
cluster. Section 4 is devoted to the derivation of the halo 
MF and §5 to the conclusion. 

2. Luminosity Functions 

Until recently, derivations of stellar LFs in GCs down to 
faint absolute magnitudes have been based on observa- 
tions made from the ground, thus suffering from crowd- 
ing effects and uncertain completeness corrections (see e.g. 
Richer et al. 1991). Data from the HST have now become 
available for several GCs down to <~ 4 magnitudes fainter 
limits (My <~ 14) without the need for large completeness 
corrections, spaning a large range of metallicity from so- 
lar to 1% solar. The MS is sharp and well defined down 
to almost the limit of the observations, which corresponds 
to m <~ 0.1 M Q (BCAH), and the LFs are essentially com- 
plete (within ~ 80%) down to this limit. These observa- 
tions are presented and discussed in §2 of BCAH, with 
the different parameters, distance moduli and reddening 
corrections, used for these clusters. 

Although there is a general agreement between the 
ground-based and the HST LFs at bright and intermediate 
magnitudes, they differ substantially at and beyond the 
magnitude limit of ground-based observations. Although 
the ground-based LFs were generally rising down to the 
limit of the observations (Mj ~ 9), most of the LFs de- 
rived with the HST show a clear decline beyond this limit. 

A comparison between the LFs of some of the afore- 
mentioned clusters, namely 47Tuc, NGC6397 and coCen, 
has been done by Cool et al. (1996) and Santiago et al. 
(1996). In spite of the differences between these clusters 
arising from their age, metallicity, galactic position, and 
thus dynamical evolution, the LFs are found to be very 
similar, down to the limit of the observations, except for 
the faint end oitoCen. 

The clusters have been observed in different passbands 
characteristic of the HST, namely F814 ~ I c , F555 ~ Vj 
and F606 where Ic and Vj denote the standard Johnson- 
Cousins (JC) system. HST colors were transformed into 
the JC system, using the Holtzman et al. (1995) synthetic 
transformations. As shown in BCAH, agreement between 
the theoretical and the observed CMDs is remarkable in 
all filters, HST and JC, assessing both the accuracy of the 
models and of the Holtzman photometric transformations 
(see BCAH for details). In Fig. 1, we show the bolometric 
LFs, corrected for incompleteness, for the afore-mentioned 
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Table I. Metallicity, distance moduli (corrected from extinctions; see caption of Fig. 1) and Galactic location of the globular 
clusters. [Fe/H] is the observed metallicity, [M/H] is the metallicity used in the models, which takes into account the enrichment 
in a-elements (see BCAH for details). Zgp, Rgc and r/ru denote the height above or below the Galactic disk, the distance 
from the Galactic center and the location of the observations w.r.t. the half-mass radius, respectively. 
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clusters, using the published LFs and the bolometric cor- 
rections of Chabrier & Baraffe (1997a) and BCAH. We 
have used cubic spline interpolations to interpolate the re- 
lations M^(Mboi) and their derivatives dMA/dMboi (global 
polynomial fits are not reliable). We verified that LFs of 
the same cluster observed in different passbands yield the 
same bolometric LF. This adds credibility to the present 
bolometric corrections and stellar models. 

We believe these results to illustrate consistent com- 
parisons between the LFs of various clusters, over a large 
range of metallicities, from near solar to [M/H] = —2.0. 
The LFs have been normalized arbitrarily to M^ Q \ — 8, i.e. 
m <~ 0.4 — 0.5 M Q over the present range of metallicity. 
For all the considered GCs, both the degree of complete- 
ness and of photometric accuracy are maximum in this 
region and stellar evolution effects do not affect this mass 
range, even for the older clusters. In Table I, we give the 
galactocentric position Rgc, the height above or below 
the Galactic plane Zqp and the metallicity for the var- 
ious clusters. These values have been averaged over the 
slightly different values found in the literature. We also 
give the location of the observations w.r.t. the half-mass 
radius r^. Some indetermination (~ 0.2 mag) remains for 
the distance modulus of various clusters, as quoted by the 
different observers. We have carried out calculations with 
the different values and found that the resulting uncer- 
tainty in the MFs remains within the observational error 
bars. The extinction corrections are the ones quoted in 
Table I of BCAH and in the different references corre- 
sponding to the observations. 



As seen in Fig. 1, all LFs exhibit a similar quali- 
tative behaviour in the region Mboi ^ 9, i.e. Mgu ~ 
Mi <, 8.5, M 555 w M v <, 10, with a rising slope 
Alog(dN/dM bol )/AM hol ~ 0.25, i.e. dN/d\og(L/L Q ) cx 
-(L/L Q )-°- 6 . Below this limit, some differences appear 
between the various clusters. Whereas some LFs, e.g. for 
coCen, exhibit a flatfish behaviour, most of them drop 
more or less rapidly at the faint end. As will be shown 
in the next section, these different behaviours, and the 
respective magnitudes corresponding to the peak in the 
LFs of the clusters, reflect the dependence of the mass- 
magnitude relationship on metallicity, from a rather uni- 
versal mass-function. The effect of metallicity in the mass- 
magnitude relationship translates into a slightly differ- 
ent location for the peak of the luminosity function, cor- 
responding to m ~ 0.2 Mq, lower metallicities yielding 
brighter magnitudes for the same mass, an intrinsic prop- 
erty of stellar matter (see BCAH). The case of iVGC6397 
deserves particular attention, for there is a substantial 
qualitative difference in the faint end of the LF for this 
cluster between the observations of Paresce et al. (1995b) 
and Cool et al. (1996). Whereas Paresce et al. get a 
strongly decreasing LF, Cool et al. argue for a milder de- 
crease, or even for a flatfish behaviour at the faint end 
of the LF. Since both observations were completed in a 
similar region of the cluster (r/r/, ~ 3), the difference can 
not stem from different mass segregation effects. A natural 
explanation would be incompleteness at faint magnitudes 
in the Paresce et al. sample but this hypothesis seems to 
be excluded by these authors. A third explanation would 
be the bias in the observed LF due to bining at faint mag- 
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Fig. 1. Bolometric luminosity functions of different globu- 
lar clusters, obtained from the observed LFs in various pass- 
bands, including F814, F555 and _F606, using the bolomet- 
ric corrections of Baraffe et al. (1997; BCAH) and Chabrier & 
Baraffe (1997a). The clusters are displayed by increasing metal- 
licity (see Table I) from top to bottom, namely : iVGC7099, 
NGC7078, 7VGC6341 from Cool et al. (1996) observations in 
the V (solid line) and Jgi4 (dashed line) bands and from De- 
Marchi & Paresce (1995b) in J 8 i4 (dotted line); 7VGC6656 from 
DeMarchi & Paresce (1997) in 7 8 i 4 ; iVGG6397 from Cool et al. 
(1996) and from Paresce et al., (1995b), same caption as above; 
wCen from Elson et al. (1995) (dot-dash line); JVGG6752 
(dash) and M3 (solid) from Fusi Pecci et al. (1996); AlTuc 
from De Marchi & Paresce (1995a) (solid line) and from Santi- 
ago et al. (1996) (dot-dash line). The distance moduli are given 
in Table I, corrected from extinctions with the corrections given 
in Table I of BCAH or quoted by the respective observers. All 
observations are normalized arbitrarily at Mboi = 8 for each 
cluster. The error bars (see the afore-mentioned references) are 
not shown for a sake of clarity. 

nitudes. This point will be discussed in detail in §4. In any 
case, this uncertainty prevents an accurate determination 
of the mass function (MF) at the bottom of the MS for 
NGC6397, and stresses the need for further observations 
of this cluster. 

These various LFs were obtained at similar location in 
the clusters, between ~ 1 and 3 times the half-mass radius 
(see Table I), assuming that the present h&li-light radius 
is close to the half-mass radius, a reasonable assumption 
(Djorgovski, 1993). In any case the region where the LFs 
have been determined is well inside the tidal radius. As dis- 



cussed below, this is the best location to obtain a PDMF 
as similar as possible to the initial MF (IMF), although 
dynamical evolution may have shaped it to some extent. 

In summary, the different LFs seem to be marginally 
similar down to Mboi ~ 10, in spite of the very different 
metallicities. Whether this apparent non-dependence of 
the LF upon metallicity will translate into a similar uncor- 
relation in the mass-function depends on the degree of cor- 
relation of the mass-luminosity relationship with metallic- 
ity, as examined below. 

3. Globular cluster mass-functions 

As already mentioned, the agreement between the ob- 
served GC main sequences and the theoretical BCAH ones 
over the entire sequences is excellent, less than 0.1 mag, 
for all clusters. In particular the models reproduce accu- 
rately the two changes in the slope of the color-magnitude 
relation, near ~ 0.4 — 0.5 Mq and below ~ 0.2 Mq, which 
reflect intrinsic physical properties of the atmosphere and 
the interior of VLMS (see Chabrier & Baraffe, 1997a). 
These comparisons asses the validity of the calibration of 
low-mass stars from the observed colors and magnitudes, 
as given in BCAH and in Chabrier & Baraffe (1997a). 

As demonstrated in the afore-mentioned papers, rela- 
tionships based on grey atmosphere models and T e g(r) 
relationships yield inaccurate, or at best unreliable, m — L 
and m — T e g relationships over the entire low-mass star 
domain, and most importantly near the hydrogen burning 
limit. This translates directly into unreliable MFs, thus bi- 
aising any attempt to derive robust conclusions. We stress 
also the importance of using an accurate representation of 
the m — L relation when deriving mass functions. Given 
the severe variations of the slope of this relation along 
the LMS range, a global polynomial fit, in particular, is 
a very poor and misleading representation. A polynomial 
fit does not reproduce all the points of the m — L relation 
and can lead to severely erroneous MFs, in particular near 
the bottom of the MS. 

The LFs discussed in §2 have been converted into MFs, 
for the appropriate metallicities, and are displayed on Fig. 
2. As shown in Fig. 2, all MFs are rising monotonically 
down to the limit of the observations, i.e. <~ 0.1 — 0.15 Mq, 
close to the hydrogen burning limit, except possibly for 
47Tuc where the MF flattens below ~ 0.2 M . This shows 
convincingly that the decrease of the luminosity function 
near the faint end does not reflect a lack of stars but in- 
stead the severe drop in the mass-luminosity relationship 
near the bottom of the MS and thus the scarcity of stars 
over large magnitude intervals for m ^ 0.2 Mq (see e.g. 
D'Antona & Mazzitelli, 1996; Chabrier & Baraffe, 1997a). 
We stress again that LFs observed in different filters, for 
a same cluster, yield similar MFs. The MFs shown in 
Fig. 2 for the different clusters are displayed by increasing 
metallicity from top to bottom, from [Fe/H] <~ —2.2 to 
[Fe/H] <~ —0.7. We recall that, given the enrichment of 
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Fig. 2. Mass- functions (in solar mass) for the different clus- 
ters with increasing metallicity from top to bottom, a in- 
dicates the slope of the power-law dN/dm oc mT a , as 
given below in brackets and shown by the dotted line for 
each cluster. Top to bottom : 7VGC6341, 7VGC7079 and 

7VGG7099 ([M/H] 2.0; a ~ 1.5); 7VGC6656 (dash-dot) 

and iVGC6397 (solid) ([M/H] = -1.5; a ~ 0.5); wCen 

([M/H] ~ -1.3; a ~ 1.0); 7VGC6752 and M3 ([M/H] 1.0; 

a - 0.5); 47Tuc ([M/H] 0.5; a ~ 0.7). Several lines for a 

same cluster indicate different observations (see Fig. 1). 



halo and globular star clusters in a-elements, an observed 
[Fe/H] for metal-poor stars corresponds to a solar-mix 
scaled metallicity [M/H] ~ [Fe/H] + 0.35 (Ryan & Nor- 
ris, 1991; see BCAH for details). A first interesting result 
is that, except for the MF of AGG6397 derived from the 
Paresce et al. (1995b) LF, all MFs are fairly well repre- 
sented by simple power laws dN/dm cx m~ a with a > 0. 
The second striking result is the weak variation of the 
slope of the MF with metallicity, with a ~ 0.5 to 1.5. 
This dependence on metallicity is pondered to some ex- 
tent by the location of the clusters in the Galaxy. This 
influence, however, remains weak as long as the clusters 
are located far enough 1 kpc) from both the disk and 
the Galactic center, as shown for example for AGG6752 
or wCen, located about 1 kpc above the Galactic plane, 
which have a slope similar to or steeper than M3, of com- 
parable metallicity, located ~ 10 kpc above the plane. In 



the same vein, M30 is found to have a slighly steeper MF 
than AGG6341, for the same metallicity, although it lies 
closer to the Galactic center. This shows that there is no 
direct link between the MF near and the location in the 
Galaxy and that some internal dynamical process, namely 
mass segregation, may have slightly modified the IMF of 
the clusters near the half-mass radius. Although the last 
bin of AGG6656 seems to indicate a decrease at the faint 
end of the MF, this feature is likely to stem from incom- 
pleteness due to crowding near the bottom of the MS (De- 
Marchi & Paresce, 1997). Interestingly enough, however, 
both the MF of AGG6656 and AGG6397, as observed by 
DcMarchi & Paresce, seem to indicate a flattening below 
~ 0.2 — 0.3 M Q (although the error bars in NGC6656 are 
consistent with a slowly rising MF), as suggested by these 
authors. Even though both clusters have similar Galactic 
positions (see Table I) , they have quite different orbits and 
thus have experienced different dynamical histories (De- 
Marchi & Paresce, 1997). This prompted these authors to 
suggest that this behaviour is an intrinsic properties of 
GC's. We note, however, that this result differs from the 
ones obtained for all the other clusters and more impor- 
tantly from observations of the same cluster (AGG6397) 
by Cool et al. (1996) (see below), as mentioned previously. 
Solving this issue requires further observations of these 
clusters near the bottom of the MS. The case of A7Tuc 
is also of particular interest, for the MF of this relatively 
metal-rich cluster seems to flatten off below <~ 0.2 M©, 
although the lack of observations beyond this limit pre- 
vents robust conclusions. The detailed study of the MF of 
metal-rich clusters and field stars is presently under work 
(Chabrier, Mera & Baraffe, 1997). 

The similitude between the MFs of the different clus- 
ters, all determined at radial distances near the half mass 
radius, seems to confirm that the stellar population in this 
region is not very sensitive at least to external dynamical 
processes. In particular, given its large distance from the 
Galactic plane and from the Galactic center, M15 is be- 
lieved to have an almost uncontaminated MF, and is thus 
of particular importance for the understanding of the dy- 
namical evolution of GCs. It is instructive to note that the 
MF of M15 (and M3) is very similar to the one of clus- 
ters located closer to the Galactic plane, although these 
clusters are expected to have experienced quite different 
dynamical histories, in particular very different tidal strip- 
ping and disk/bulge shocking, resulting in a substantially 
different fraction of VLMS at the periphery. 

As mentioned in the previous section, the indetermina- 
tion near the faint end of the LF for AGG6397 translates 
into a difference in the MF below ~ 0.3 Mq. This is il- 
lustrated in Figure 3. The two LFs of Cool et al. (1996) 
(in the F555 and F814 filters) yield a monotonic, rising 
MF, whereas the LF of Paresce et al. (1995) yields a flat- 
fish MF below 0.3 Mq, although with large error bars 
near the faint end. The difference remains when using 
any of the two distance moduli quoted for this cluster 
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Fig. 3. Mass function for 7VGC6397 from different observed 
LFs : solid line : Jsi4 (Cool et al., 1996); dashed line : V555 (Cool 
et al., 1996); dot-dashed line : Jsi4 (Paresce et al., 1995). Also 
shown is the MF obtained from the LF observed by Mould 
et al. (1996) at r ~ 6 x r h (dotted line) to illustrate mass 
segregation effects. The MFs are normalized at 0.5 Mq. 

and stems from the very difference at the faint end of the 
two LFs. For this reason, we will leave this cluster aside 
for the rest of the discussion until more observations can 
nail down this issue. Note that the low galactic position 
of NGC6397, Z w 0.5 kpc, its proximity to the Galac- 
tic center, and its likeky collapsed state suggest that this 
cluster has experienced severe disk and bulge shocking, 
resulting in an accumulation of low-mass stars at the pe- 
riphery and a substantial depletion even near rh- Note 
also that the relatively short half-mass relaxation time of 
NGC6397 (t h ~ 3 x 10 8 yr; Djorgovski et al., 1993) may 
have led already to substantial evaporation in this region. 
This might provide a simple explanation for the MF in- 
ferred from Paresce et al. data and to explain the fact that 
NGC6397 has a shallower MF, and thus a smaller amount 
of LMS than clusters of similar metallicity, e.g. M30. On 
the opposite, if the MF keeps rising near the bottom of 
the MS, as inferred from the Cool et al. LF, the similarity 
of the MF with the ones of clusters located far from the 
Galactic plane shows that even strong tidal stripping and 
disk shocking does not strongly affect the inner regions of 
the cluster, near rh- 

An interesting information about tidal stripping and 
mass segregation in NGC6397 arises from the observa- 
tions of Mould et al. (1996) who obtained the LF of the 
cluster with the HST at a larger radial distance than the 
aforc-mcntioncd observations, namely at 10 arcmin from 
the center, i.e. r/r^ <~ 6. The MF obtained from Mould et 
al. LF is shown in Fig. 3 (dotted line). The steeper slope 



(a <~ 1.3 instead of a ~ 0.5) and thus the larger number 
of very low-mass stars (SS 0.5 Mq) w.r.t. to the previous 
determinations closer to r/, illustrates the onset of tidal 
striping and mass segregation, yielding an excess of very- 
low-mass stars in the outer parts of the cluster. Therefore 
it is not excluded that VLMS near rh may already have 
been stripped in this cluster. In any case, it shows that no 
robust conclusion about the IMF of GCs can be inferred at 
present from the observations of NGC6397, and probably 
not from any cluster located at low Galactic latitude. 

The core relaxation times of the present clusters are 
~ 10 8 — 10 9 yr, at least one order of magnitude less that 
the ages of the clusters, so that we expect the stars near 
the central core of the clusters to be close to thermal equi- 
librium. Therefore low-mass stars in the core should be 
substantially depleted as a result of mass segregation. The 
half-mass radius is in general about 50 times larger than 
the core radius, so that we expect mass segregation to 
be mild in this region. The effect of mass segregation in 
GCs has been demonstrated by various investigations of 
the core (r c w 0.01 rh) of different clusters (DeMarchi & 
Paresce, 1996; Paresce et al., 1995a; King et al., 1995). The 
severe deficiency of faint stars in the center of the clusters 
translates into a strongly decreasing MF below ~ 0.6M Q . 
This reflects the effect of mass segregation in the central 
core, due to two-body relaxation, as expected from nearly 
thermodynamic equilibrium, providing the time for core 
collapse is much larger than the half-mass relaxation time, 
which is the case in general. This result is confirmed by 
the analysis of King, Sosin & Cool (1995) who examined 
the relation between the local MF of the cluster, as de- 
termined presently, and the global MF, using a dynamical 
model for the evolution of the cluster. As shown by these 
authors, whereas the MF determined in the central core is 
severely affected by mass segregation, the MF determined 
near the half-mass radius is quite similar to the global MF. 

3.1. Predictions in the NICMOS filters 

The observations of GCs in near infrared colors will soon 
be possible with the HST NICMOS camera, in the F110VF 
(~ J), FimW (~ H) and F187W filters, where J and 
H are defined in the CIT system. Mass-magnitude re- 
lationships and CMDs have been derived by BCAH in 
these niters. In Figures 4a-c, we show the predicted lu- 
minosity functions in these filters for [M/H] = —2, as 
obtained with three different rising MFs (a > 0) and a 
flat MF (a = 0) beyond 0.2 M . The LFs are normal- 
ized to iVtot = *" ax (dN/dm) dm — 1 , where N to t is the 
total number of stars between 0.8 Mq, about the turn- 
off mass, and the hydrogen burning limit, which corre- 
sponds to magnitudes M m ; n and M max , respectively, for 
the appropriate metallicity, as given in BCAH. The slope 
of the MF, at least for a > 0.5, could be determined by 
future observations from differential comparisons of the 
star counts on each side of the peak of the LF, provid- 
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Fig. 4. Predicted LFs for clusters with [M/H] = — 2 for different mass functions : dotted line : a = below 0.2 Mq; dashed 
line : a = 0.5; solid line : a = 1.0; dot-dash line : a = 1.5 in the NICMOS filters F160W (Fig. 4a), F187W (Fig. 4b) and F110W 
(Fig. 4c). The effect of metallicity is illustrated on Fig. 4d for a = 1.5 : [M/H] = -0.5 (dotted line), [M/H] = -1.0 (dot-dash 
line), [M/H] = -1.5 (solid line) and [M/H] = -2.0 (dashed line). 



ing accurate completeness corrections. Steeper MFs yield 
more (resp. less) VLMS below (resp. above) ~ 0.2 Mq, 
which corresponds to the bright peak in the LFs. A flat 
MF yield a larger number of stars and a small plateau 
between <~ 0.1 — 0.2 Mq, followed by an abrupt drop. The 
effect of metallicity on the LF is illustrated in Figure 4d 
for a — 1.5. Although the qualitative behaviour of the 
LF remains the same, the effect of metallicity is a global 
shift of the LF toward slightly brighter magnitudes with 
decreasing metallicity, a direct consequence of the metal- 
licity dependence of the mass-magnitude relationship (see 
BCAH, Chabrier & Baraffe, 1997a). The effect, however, 



remains small and shows the weak dependence of the LF 
upon metallicity, at least for metal-depleted objects. 

An other interesting feature about these LFs is the ex- 
pected brown dwarf signature, as illustrated by the peak 
at faint magnitudes in Figs. 4a-c. The gap in the LF stems 
from the large luminosity drop in the mass-luminosity re- 
lationship in the stellar to substellar transition domain, 
which corresponds to several magnitudes over a few hun- 
dreths of a solar mass, whereas the sharp increase reflects 
the flatter slope of this relation in the brown dwarf regime. 
We stopped our calculations at m — 0.055 Mq, which 
corresponds to the coolest non-grey model atmospheres 
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presently available. This is enough to illustrate the un- 
missing signature of even a moderate BD population in 
GCs, as obtained with the mentioned MFs. Interestingly 
enough, we note that the gap in the LF varies significantly 
with the filters, which reflects the strong variation of the 
BD energy distribution in the near infrared frequency do- 
main. Even in the most favorable case, i.e. F110, the gen- 
uine brown dwarf domain lies ~ 4 mag below the bottom 
of the main sequence. This sets the scale for BD search in 
GCs. 

3. 2. Brown dwarf mass fraction 

The mass under the form of brown dwarfs in GCs is de- 
termined by : 

m d^ dm = 2^ (W 'P Q " m *™ ] (1) 

where A is the normalization of the MF dN/dm = Am~ a . 
We have normalized all the MFs at ~ 0.4 — 0.5 M Q , as 
mentioned previously (see §2). The value m up = 0.083 M 
is the hydrogen burning minimum mass for [M/H] — — 2 
and about it for metal-poor stars in general (Chabrier & 
Baraffe, 1997a) whereas m\ ow is the minimum mass of 
the brown dwarf IMF. Since this quantity is presently 
unknown, the maximum brown dwarf contribution cor- 
responds to m\ ow = 0, since the integral is convergent 
(a < 2). A more realistic value mi ow = 0.01 M@, which 
corresponds to the deuterium burning minimum mass 
(Saumon et al., 1996) and to the minimum Jeans mass at 
zero-metallicity (Silk, 1977) does not change the results 
significantly, given the small contribution of the lowest 
masses for weakly rising MFs. 

Even the steepest MF a = 1.5 yields a brown dwarf 
contribution ~ 30% of the stellar mass fraction below 0.8 
Mq, about the cluster turn-off mass. A flat MF or a larger 
minimum mass m\ ow yield essentially no mass (< 1%) in 
the brown dwarf domain. These calculations show that, 
even though there might be a substantial number of BDs 
in GCs (rising MF), their mass fraction is fairly small. 
This confirms that there is essentially no dark matter in 
GCs. 

4. Halo mass function 

If the PDMF slope is not strongly correlated with the 
metallicity, we expect the IMF for GCs to be rather similar 
to that of the disk and the halo. This issue is addressed 
below. 

As shown in §4.3 of BCAH, there is no apparent dif- 
ference in ±he sequences of globular star clusters and halo 
field starsCL This corroborates the traditional view that 

* The term halo is used for different populations in the litera- 
ture and leads to some confusion. In the present paper, we de- 



globular clusters and halo field stars are members of the 
same stellar population (see e.g. Fall & Rees, 1985). 

The direct determination of the halo LF is a tremen- 
dously difficult task. A 0.09 M© star with [M/H] = -0.5, 
about the metallicity of the old (thick) disk population 
(see e.g. Leggett, 1992), has a magnitude My ~ 17 
(Baraffe et al., 1995). Since the old disk scale height is 
~1-1.5 kpc (Kuijken & Gilmore, 1991; Bienayme, Robin 
& Creze, 1987), we expect the disk population to con- 
tribute significantly to star counts up to at least V ~ 27, 
I ~ 23. This clearly illustrates the difficulty of an ac- 
curate determination of the spheroid or halo LF from 
ground-based photometric measurements. This difficulty 
is in principle circumvented with the Hubble Space Tele- 
scope, which reaches V ~ 30 in the Hubble Deep Field 
(HDF) with 0.1 arcsec resolution, but the small field of 
the HDF (~ 5 arcmin 2 ) yields too small statistics to de- 
rive a reliable halo LF. 

This is why we must rely on the nearby halo LF, based 
on geometric parallax determinations of high- velocity faint 
stars in the solar neighborhood, as obtained by Dahn et 
al. (1995). The high tangential velocity of the stars in 
the Dahn et al. LF («tan > 220 km.s -1 ) is a strong in- 
dication of a halo population. As shown in BCAH (Fig. 
8), the sequence of high- velocity subdwarfs of Monet et 
al. (1992), part of the Dahn et al. (1995) sample, is fully 
consistent with the GC sequences with an average metal- 
licity [M/H] ~ -1.5 and a dispersion A[M/H] m ±0.5 
(BCAH). However, as shown in the same afore-mentioned 
figure of BCAH, a substantial fraction of the subdwarfs 
observed by Dahn et al. (1995), which all have halo- 
like kinematic properties, exhibit colors characteristic of 
metal-rich ([M/H] *S —0.5) objects. The very nature of 
these stars remains a puzzle. Given the extremely low 
probability for a large number of disk objects to have such 
large velocity dispersions, these objects are likely to be 
halo stars formed in the Galactic central region. A con- 
sistent, bin by bin analysis of these observations requires 
precise completeness corrections, which are not available, 
so we can only derive a global analysis. The luminosity 
of LMS increases with decreasing metallicity, for a given 
mass (BCAH). Conversely, for a given magnitude, in the 
range of interest, the mass of a LMS with [M/H] = —0.5 
is —0.1 M Q larger than with [M/H] = —1.5. Therefore, 
assuming a metallicity [M/H] = — 1.5 for all the stars en- 
tering the Dahn et al. LF will slightly overestimate the 
low-mass end of the MF and thus provides an upper limit 
for the slope and the normalization of the halo MF at the 
bottom of the MS. 



note stellar halo or spheroid the stellar population determined 
by the Hubble (l/r A ) or the DeVaucouleurs (e~ r ) density 
profile, whereas we call dark halo the outer 1/r 2 density distri- 
bution. We will use the generic term halo to denote the overall 
population belonging to the afore-mentioned spheroid and dark 
halo, in opposition to the disk population. 
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The MF determined from the Dahn ct al. LF for this 
metallicity is shown in Fig. 5 (solid line). Also shown is the 
MF derived from the photometric spheroid LF determined 
by Richer & Fahlman (1992, RF) (dashed line). As shown 
in Fig. 5 (see also Mera, Chabrier & Schaeffer, 1996), both 
MFs arc fairly consistent for m ^ 0.15 M Q but differ sig- 
nificantly below this limit. The RF's MF rises abruptedly 
whereas the Dahn et al. MF drops below m <~ 0.10 M©. 
Such a qualitative difference is surprising, for both LFs are 
supposed to probe the same stellar population, essentially 
the spheroid one. As mentioned above (see also Reid et al., 
1996), the last bin in the RF LF is likely to be contami- 
nated by old disk stars and can not be used to derive a halo 
MF. The difference between the two MFs shown in Figure 

5 becomes significant below m ~ 0.15 M©, i.e. My <; 12, 
Mi £ 10, M bol £ 11 for [M/H] = -1.5 (BCAH; Chabrier 

6 Baraffe, 1997a). At 15 kpc this corresponds to V ~ 28, 
I ~ 26. At faint magnitudes, e.g. V ~ 25, an old-disk 
0.4 M© star {[M/H] > -1.0; M v ~ 10) at 10 kpc can 
be misidentified as a metal-poor <~ 0.1 M© halo star 
{[M/H] < -1.5; My - 15) at 1 kpc with similar V -I. 
This will yield a severely overestimated halo number den- 
sity near the bottom of the MS. Note also that there are 
only 3 stars in the two last bins of RF's LF. 

On the other hand, the last 3 bins in the Dahn et al. 
LF rely on a fairly small number of stars (7, 5 and 1, re- 
spectively), and have a low statistical significance. At last 
we want to emphasize the bias introduced by bining the 
data at faint magnitudes. The number of stars per unit 
magnitude is about 5 times smaller below ~ 0.1 M© than 
above this value (see e.g. BCAH), because of the abrupt 
drop in the mass-luminosity relationship. A good statis- 
tics thus implies a large magnitude interval but this in turn 
yields a large indetermination in the exact magnitude of 
the stars belonging to the interval. The faint end of the 
Dahn et al. LF is likely to be hampered by these limita- 
tions, due mainly to incompleteness and to the scarcity of 
objects in a large magnitude range at the very bottom of 
the main sequence. Note that the same observational bias 
could apply to the uncertain determination of the faint 
end of 7VGC6397 LF, as discussed in the previous section. 

As shown in Fig. 5, the Dahn et al. MF is consis- 
tent with a rising MF dN/dm oc m~ a with a <~ 2 down 
to 0.12 M© at the 1-a (99%) confidence level, and with 
a <~ 1.6 down to the very last bin, i.e. 0.09 M© at the 
2-cr (90%) C.L., for [M/H] = -1.5. This C.L. increases 
if we use a slightly lower metallicity. A x 2 -test applied to 
the MF determined in Figure 5 yields a most-likely slope 
a ~ 1.64, assuming a power-law MF, although we cannot 
exclude a non power-law behaviour at the very end of the 
MF. Interestingly enough, similar conclusions have been 
reached recently by Richer & Fahlman (1996) who reana- 
lyzed their photometric LF and derived a MF with a slope 
a ~ 1.7. A more precise determination of the shape of the 
MF near the hydrogen-burning limit, which bears essen- 
tial consequences for star formation theory, requires larger 
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Fig. 5. Mass function for halo field stars. Solid line : MF deter- 
mined from the nearby LF of high-velocity subdwarf (Dahn et 
al., 1995); dashed line : MF determined from the photometric 
LF of Richer & Fahlman (1992), both for [M/H] = -1.5. Dot 
and dash-dot lines : power-law MFs with a = 1.6 and a = 2.0, 
respectively. 

statistics in this region. However the present determina- 
tion is sufficient to determine with reasonable accuracy the 
halo stellar and substellar mass fractions, as examined in 
the next section. 

5. Halo stellar and substellar mass fractions 

The spheroid extends between ~ 10 and 20 kpc from the 
Galactic center (see e.g. Reid ct al., 1996 and below), i.e. 
V ~ 25-30, 1 ~ 23-28 for the bottom of the MS of metal- 
poor stars (BCAH). This is out of reach with ground- 
based observations but can be probed partially with the 
HDF, which reaches V ~ 30, I ~ 28 (Flynn, Gould & 
Bahcall, 1996). As shown in the previous section, a power- 
law MF dN/dm oc mT a with a slope a ~ 1.6 — 2.0 is fairly 
consistent with all halo MFs (excluding the last bin of RF) 
down to at least ~ 0.12 M©, at the 99% C.L., and down to 
the limit of the observations at the 90% C.L. This yields 
the following halo MF : 

/x fc (m) « 4.0 ±2 x io-*(-JZL-)-(i-™*) M-V- 3 (2) 

Integration of this MF gives the halo stellar (m < 
0.8 M© for t > 10 Gyr) density in the solar neighborhood: 

PK « 1.0 x 10" 4 M .pc" 3 (3) 

Comparison of this value with the disk low-mass star 
density (Mera, Chabrier & Baraffe, 1996), and assuming 
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that the MF (2) is valid for masses up to m = 0.8 Mq 
(My ~ 4) gives the mass- normalization of the halo/disk 
total stellar population (including the contribution of 
more massive stars for the disk population) down to the 
bottom of the MS, p h J Pd* ~ 1/400. 

5.1. Star count analysis 

Although the small field of the Hubble Deep Field prevents 
the determination of an accurate halo LF, as mentioned 
previously, it provides a stringent constraint on the halo 
mass function and mass model. 

The number of stars per arcmin 2 between apparent 
magnitudes e.g. I m i n and I max toward Galactic coordi- 
nates (b, I) reads : 



dN_ 
~dtt 



3600 v 180 ; 



olM, 



dr r 2 p(r(l, b)) 



alN dm 
dm dMj 



(4) 



where (i max /10pc 



10" 



(resp. d m i n ) denotes the 



maximum (rcsp. minimum) distance to which a star of 
magnitude Mj will be visible, M m i n -M max the magnitudes 
corresponding to the mass range of interest (< 0.8 Mq) 
and p(r) is the mass density profile characteristic of the 
Galactic region considered. The mass-magnitude relation 
dm/dMj is taken from BCAH. Table II gives the pre- 
dicted star counts (between 0.8 Mq and the hydrogen 
burning limit) at faint magnitudes in the Hubble Deep 
field, I = 126°, b = 55°, for the observed field of view, 
0.00157 sq. degrees. The Galactic model consists of a thick 
disk with an exponential oc e~l z l/ z ° or a a sech 2 -^ pro- 
file with different scale heigths, zq = 1000 and 1500 pc, 
a spheroid with a DeVaucouleurs oc e - w ( R / R <sY 4 profile 
and a flattened (q = 0.6) dark halo with an isothermal 
distribution p(r) oc 1/r 2 . The mass functions are the ones 
determined in Mera et al. (1996b) for the disk and in Eq. 
(2) for the halo. The average metallicity is [M/H] = —0.5 
for the disk and [M/H] = —1.5 for the halo, respectively. 
As shown in Table 2, the results obtained with spheri- 
cally symetric distributions are in severe conflict with the 
observations, whereas the observed star counts are consis- 
tent with an oblated q — c/a <~ 0.6 spheroid plus a sech 2 
thick disk with a scale height zq & 1 kpc. The observations 
of the HDF exclude completely the presence of a signif- 
icant dark halo stellar population. Therefore, almost all 
the stars in the Dahn et al. sample belong to the spheroid. 

An event more stringent comparison is given in Ta- 
ble III which compares the predicted star counts with 
the HDF observations in two different passbands, namely 
the V-band (Mendez et al., 1996) and the I-band (Reid 
et al., 1996). The Galactic model consists of the afore- 
mentioned thick disk characterized by a density profile 



p(r) = 0.05 pQe-( R - R °V l ° d sech 2 ^ with a local den- 
sity pq s=s O.IMq pc" 3 , a solar galactocentric distance 
Rq — 8.5 kpc, a scale length l d = 2.8 kpc and a scale 
height zq = 1 kpc, and a flattened (q = 0.6) DeVau- 
couleurs spheroid. The spheroid and disk populations are 
determined from their location in the CMD (see e.g. Reid 
et al., 1996 Fig. 13; Baraffe et al., 1995). We note the 
excellent agreement (within ~ lcr) between theory and 
observations in the two passbands, for every color and 
magnitude interval. The spheroid VLMS population ex- 
tends up to I <~ 27, V ~ 29 and vanishes beyond this 
limit. Most of the spheroid stars in the range 25 < I < 27, 
1.5 ;S> V — I ^ 2.5 have masses between 0.1 and 0.2 Mq, as 
determined from the present stellar models and mass-color 
relations for different metallicities (BCAH). The scarcity 
of stars at fainter magnitudes thus reflects the bottom of 
the spheroid MS. The absolute magnitude of a 0.1 M 

star with [M/H] 1.5, -2.0 is My ~ 13.5, Mj - 11.3 

(BCAH) so that these stars are located at ~ 10 — 15 kpc 
from the Sun. The present analysis provides a relatively 
precise determination of the size of the Galactic spheroid. 
It also shows the consistency of the present disk and halo 
mass-functions and Galactic model with HDF faint mag- 
nitude observations in different colors. It also confirms 
the negligible stellar population in the dark halo, at most 
~ 1% of the spheroid population, as inferred from Table II. 
The present analysis shows that we have seen essentially 
all main sequence stars in the Galaxy. 

The spheroid stellar density (3) is thus about 1% of the 
dynamically-determined dark matter density and 0.1% of 
the disk total density, whereas the dark halo stellar density 
is about two orders of magnitude smaller, as implied by 
the present star-count analysis. This yields a contribution 
of the spheroid to the optical depth towards the LMC 
T sp h ~5x 10~ 9 , the eventual dark-halo contribution being 
negligible. 

The extrapolation in the BD domain, as done in the 
previous section, yields a negligible amount of substellar 
objects in the Galactic halo. A minimum mass m\ ow — 
0.01 Mq for the BDs and the extrapolation of the MF (2) 
yields a BD number-density in the spheroid: 



n hBD ~ 2± 1 x 10 d pc ' 
and a BD mass-density : 



Phi 



0.5 ±0.5 x lO~ 4 M pc~ 



(5) 



(6) 



i.e. ^ 1% of the dark matter density. This is consistent 
with a detailed analysis of the MACHO+EROS microlens- 
ing observations toward the LMC (Alcock et al, 1996; 
Aubourg et al., 1993; Mera, Chabrier & Schacffer, 1997). A 
substantial fraction of halo brown dwarfs not only would 
imply a MF below the hydrogen burning limit radically 
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Table II. Different contributions to the star counts in the 5.6 Arcmin 2 Hubble Deep Field (b,l) = (55°, 126°) (Reid et al., 
1996) with the mass function determined in Mera et al. (1996b) for the disk and in Eq. (2) for the halo. We considered different 
density-profiles (see text) .The repartition between disk and halo in the observations is just indicative and has been obtained 
from color criteria (see e.g. Reid et al., 1996 Fig. 13). 



I = 20 - 26 / = 24 - 26 



old disk 



e- |z|/zo ; 2 = 1.0 kpc 
e- |z|/zo ; 2 = 1.5 kpc 
sech 2 (^-); z = 1.0 kpc 
sech 2 (^); z = 1.5 kpc 


9 

23 
3.9 
12.2 


2.6 
8.7 
0.7 
3.2 


Observations 


4-6 


1-2 


spheroid 

q = l 
q = 0.6 


12.3 
5.8 


5.4 
1.8 


dark halo 

q = 0.6 


244 


176 


Observations 


5-7 


3-4 



Table III. Star counts in the Hubble Deep Field in the V-band (Mendez et al., 1996; 4.69 arcmin 2 ) and in the I-band (Reid 
et al., 1996; 5.6 arcmin 2 ) (tabulated CMDs). See text for the Galactic model. The average metallicity is [M/H] = —0.5 for the 
disk and [M/H] = —1.5 for the spheroid, respectively. The numbers in brackets for the spheroid correspond to [M/H] — —2. 
At fainter magnitudes, the predicted star counts are consistent with zero in the color-range of interest. 



i / (v-i) 



0.5-1.5 



1.5-2.5 



2.5-3.5 



18-20 

Thick disk 0.1 0.8 0.2 

Spheroid 0.4(0.5) 0.1(0) 0(0) 

Obs. 1 2 
20-22 

Thick disk 1.4 0.3 

Spheroid 0.6(1) 0.6(0.3) 0(0) 

Obs. 2 1 1 
22-24 

Thick disk 0.6 0.9 

Spheroid 0.3(0.8) 1.7(1.3) 0.1(0) 

Obs. 2 2 1 
24-26 

Thick disk 0.3 

Spheroid 0(0.2) 1.6(1.5) 0.3(0.2) 

Obs. 2 1 



V / (V-I) 0.5-1.5 1.5-2.5 2.5-3.5 



21-23 



23-25 



25-27 



1.2 0.1 

0.5(1) 0.5(0.3) 0(0) 

2 

0.9 0.6 
0.3(0.8) 1.5(1.1) 0(0) 

1 2 1 

0.1 0.7 
0(0.2) 1.6(1.4) 0.1(0.1) 

1 1 



different (much steeper) from the one determined by Eq. 
(2), but would also imply a large rate of short-time (SS 20 
days) events (r oc 1/y/m), whereas no event beyond 19 
days has been observed by either the MACHO or the 
EROS experiments towards the LMC. It confirms that 
the last bin of the RF LF, which would yield a ~ 5 is not 
correct. This clearly excludes brown dwarfs as plausible 
candidates for the halo missing mass and for the observed 
events towards the LMC. 



Interestingly enough, we note that the MF (2) has 
a slope very similar to the one determined previously 
for the disk population down to the bottom of the MS, 
1.5 ^ a ^ 2.5 (Mera, Chabrier & Baraffe, 1996). The 
small number statistics at faint magnitudes prevents a 
more precise determination. Note that such a rising MF 
is consistent with microlensing observations towards the 
central bulge whereas a decreasing MF near the hydrogen 
burning limit would be in severe conflict (> 5er) with these 
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observations (Han & Gould, 1996; Mera et al., 1997). If 
confirmed by future observations, this strongly suggests 
that the stellar MF is a rather universal property, almost 
independent on the metallicity, and reflects the nature of 
a dominant physical process in star formation. It would 
also show that a Salpeter-like MF extends to about three 
orders of magnitude in mass, then justifying the extrapo- 
lation into the brown dwarf domain. Note that the metal- 
licity might be detcrmincnt for the minimum mass of the 
MF. 

6. Conclusion 

In this paper, we have first transformed the LFs of several 
globular clusters observed with the HST in various pass- 
bands down to almost the bottom of the main sequence 
into bolometric LFs. All these LFs exhibit a similar be- 
haviour with a peak around Mboi ~9-10 and a drop be- 
yond. This reflects the sharp drop in the mass-luminosity 
relationship below <~ 0.2 Mq. We have converted these LFs 
into mass functions. We emphasize that these mass func- 
tions rely on mass-luminosity relationships obtained with 
stellar models whose accuracy has been demonstrated un- 
ambigously (BCAH). Since the luminosity functions for 
these clusters have been determined near the half-mass 
radius, the present-day local mass functions arc likely to 
reflect reasonnably well the initial global mass functions. 
All the MFs are rising down to <~ 0.1 M , with the ex- 
ception of the one based on the observations of NGCQ397 
by Paresce et al. (1995) and possibly the one for 47Tuc 
which seem to flatten off below ~ 0.2 Mq. This stresses the 
need for further observations of these clusters to determine 
whether the effect is real or whether it stems from vari- 
ous observational bias at faint magnitudes, as suggested, 
for NGC6397, by comparison with the observations of the 
same cluster by Cool et al. (1996). All these MFs are well 
described by a power-law dM/dm cx mT a , with essen- 
tially no dependence of a on metallicity or the Galactic 
location, with a varying between a <~ 0.5 and a ~ 1.5. 
The clusters thus seem to exhibit the same universal ini- 
tial mass function, shaped to some extent (within about 
one unit for the slope a) by further dynamical evolution. 
The effect of metallicity appears in the mass-magnitude 
relationship and translates into a glogal shift of the LF, 
and thus of the location of the peak, which corresponds 
tom~ 0.2 M , to slightly brighter magnitudes for lower 
metallicity. Given the completely different dynamical his- 
tory of these clusters, this confirms that the stellar pop- 
ulation near the half-mass radius is weakly affected by 
cither external or internal dynamical effects. Such a con- 
clusion would have to be reconsidered if the observations 
of Paresce et al. (1995) for NGC6397 are confirmed. In 
that case, the Galactic latitude and the Galactocentric 
distance of GCs might play a determinent role below <~ 1 
kpc. These weakly rising MFs suggest that the globular 
clusters contain a substantial number of brown dwarfs but 



a negligible amount of mass, and thus of dark matter, un- 
der the form of these objects. 

The mass segregation effect is illustrated by the severe 
depletion of very low-mass stars in the central core of these 
clusters, as reflected by the dropping MF below ~ 0.6 Mq, 
and by the excess of these objects, as shown by the sub- 
stantially steeper MF, in the outer field of NGC6397. 

We have derived predicting LFs in the NICMOS fil- 
ters, corresponding to either slightly rising or flat MFs for 
the different clusters, down to the brown dwarf domain. 
Comparison between these predictions and near-future ob- 
servations should determine accurately the MF of globular 
clusters all the way down to the hydrogen burning limit 
and the brown dwarf population. 

We have applied these calculations to the determina- 
tion of the halo mass function, from the observed luminos- 
ity functions of high- velocity low-mass stars in the stellar 
neighborhood and from the direct photometric determi- 
nation at faint magnitudes. Both MFs agree fairly well 
down to ~ 0.15 M Q . Such a consistency adds support to 
the fact that they probe the same stellar population. Both 
type of observations lead to the same rising MF down to 
~ 0.15 Mq, but to two different behaviours below this 
limit. However, both LFs at the faint end are hampered 
by different observational bias and most likely are consis- 
tent with the same, slowly rising MF all the way down 
to the hydrogen-burning limit with a slope a <~ 1.7 ± 0.2, 
although a decreasing MF below ~ 0.15 M Q cannot be ex- 
cluded with presently available observations. The drop in 
the LF reflects the scarcity of objects over large magnitude 
intervals near the very bottom of the main sequence. 

Comparison with star counts at faint magnitudes ob- 
tained with the Hubble Deep Field are consistent with the 
afore-mentioned halo MF and with a sech 2 old disk with a 
scale height z$ Si 1 kpc, which contributes at least ~ 50% 
of the LF down to I <~ 25, V ~ 27. These comparisons 
yield an independent determination of the size <~ 10 — 15 
kpc and the oblatcncss q ~ 0.6 of the spheroid. These cal- 
culations show that the observed high- velocity subdwarfs 
belong to the Galactic spheroid, which dominates the LF 
up to y ~ 29, I ~ 27, and that there is essentially no 
main sequence star in the dark halo. Comparison between 
future deep field (7 > 28) observations and these calcula- 
tions should determine unambigously the dark halo stellar 
and substellar mass content. 

These calculations yield a consistent determination of 
the spheroid MF down to ~ 0.1 M Q , and its normalization 
w.r.t. the disk one. They also give an upper bound for the 
dark-halo normalization l%p sp h). They yield the con- 
tribution of the spheroid and the dark-halo contribution 
to the dark matter density in the solar neighborhood, less 
than ~ 1% and 0.01%, respectively. Extrapolating these 
mass functions into the brown dwarf domain reveals the 
negligible contribution of substellar objects to halo dark 
matter, which is consistent with the lack of short time 
events in the microlensing observations towards the LMC. 
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The mass functions determined for globular star clus- 
ters and for halo field stars have the same qualitative and 
semi-quantitative behaviour, which seems to confirm the 
common origin of these stellar populations. More strik- 
ingly, these MFs are similar to the one determined for the 
disk stellar population down to the bottom of the MS. 
This general behaviour of the mass function for different 
populations down to the hydrogen burning limit, in partic- 
ular the weak correlation with the metal content, seems to 
reflect a universal feature of all stellar systems and thus 
is likely to stem from a fundamental mechanism in the 
star formation process. Although the precise determina- 
tion of the shape of the MF requires more statistics, the 
present determination represents already a solid ground to 
infer the halo stellar and substellar moss content, as done 
in the present study. Note that the effect of metallicity 
might come into play in the determination of the minimum 
mass for the formation process. Although this minimum 
mass seems to be close to ~ 1 M for dark halo objects, 
as inferred from the microlensing observations (Chabrier, 
Segretain & Mera, 1996; Adams & Laughlin, 1996), which 
is consistent with the present calculations which exclude 
a significant low-mass star population in the dark halo, 
it lies well below this limit for spheroid and disk objects 
(Mera et al., 1997). 

The bolomctric LF's displayed in Fig. 1 are available 
at CDS via anonymous ftp: 

ftp cdsarc.u-strsbg.fr (130.79.128.5) 
or via http://cdsweb.u-strsbg.fr/Abstract.html 
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